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APPLICATION OF THE METHOD OF STEIN AND SANDERS TO THE 
CALCULATION OF VIBRATION CHARACTERISTICS OF 
A 45° DELTA- WING SPECIMEN 
By John M. Hedgepeth and Paul G. Waner, Jr. 

SUMMARY 


Generalized influence coefficients are calculated by the method of 
NACA TN 5^40 for a large-scale, built-up, 45° delta-wing specimen. These 
are used together with appropriate generalized masses to obtain the 
natural modes and frequencies in symmetric and antisymmetric free-free 
vibration. The resulting frequencies are compared with those obtained 
experimentally and are found to be consistently high. Possible sources 
of the disparities are discussed. 


INTRODUCTION 


The increased importance of flutter and other aeroelastic phenomena 
has made the accurate determination of structural stiffness properties 
a necessary part of aircraft design. For this reason, the Langley 
Structures Research Division has been conducting a program of experimental 
and theoretical research of the deflection and vibration properties of 
built-up wing structures. One of the test structures is the large-scale 
45° delta-wing specimen described in reference 1. Reference 1 also con- 
tains the details of the static and vibration tests. The experimental 
data obtained from this specimen have been used to assess the accuracy 
of two analytical methods: that of Levy (ref. 2) and the one proposed 

by Stein and Sanders (ref. 3 ). These assessments are summarized in ref- 
erences 4 and 5* 

The purpose of the present paper is to report the details of the 
calculations made by the Stein-Sanders method. Influence coefficients 
are first computed by the procedure outlined in reference 3* These coef- 
ficients are then used in conjunction with the mass properties to find 
the natural modes and frequencies of free-free vibration by means of 
matrix iteration. 
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SYMBOLS 


a k 


All quantities are in pound- inch-second, units. 

stiffness coefficient for covers (k = 0, 1, 2, 3, 4) 


BoO^Ol^u parameters defined by equation (15) 

CqCybc^y) trailing- and leading-edge coordinates, respectively 


E 

g 

[ 8 ] 

l 


m 

m k 

E k 

N 


P 0' P l' P 2 


Young's modulus of elasticity 
acceleration due to gravity 

generalized influence-coefficient matrix 
semispan 

mass per unit area 

generalized distributed mass (k = 0, 1, 2, 3j 4) 
generalized concentrated mass (k = 0, l, 2, 3, 4 ) 
tip-station index, i/e 

load, moment, and second moment about y-axis 


w 


x,y 

x s 

Ps 


deflection 

coordinate system (fig. l) 
x- coordinate of span 

spar or stringer stiffness parameter 


7 

£ 


rib stiffness parameter 
spacing between stations 
parameter defined by equation (23) 
potential energy of loads 
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de f le c t ion , s lope , 
structure) 

and half- curvature at y-axis 

( supported 


deflection, slope, 
free structure) 

and half -curvature at y-axis 

(free- 

O) 

circular frequency 



[] 

rectangular matrix 



1 1 

column matrix 



LJ 

row matrix 




DESCRIPTION OF SPECIMEN 


A schematic drawing of the delta-wing specimen is shown in figure 1. 

It has a semispan of 112 inches, a root chord of 9 6 inches, and a leading- 
edge sweep of 45°. The specimen is constructed entirely of 2024 aluminum 
alloy which is assumed to have a Young's modulus of 10.6 X 106 psi and 
a Poisson's ratio of 1 / 3 . The cover sheets are relatively thin and are 
stiffened by numerous stringers. There are four spanwise spars and a 
leading-edge spar. Closely spaced ribs provide chordwise stiffening. 

The dimensions shown in figure 1 are nominal and were used in the 
calculations for locating stations, spars, and ribs. Precise dimensions 
and details of the construction and of the weight distribution of the 
specimen are given in reference 1. Points of interest to be noted are: 

the depth tapers in the spanwise direction from inches at the carry- 

3 ^ 

through section to 1 — inches at the tip but is constant in the chordwise 

direction. With the exception of the leading-edge spar, all spars are 
unspliced for their entire length and are reinforced with heavy caps. 

The leading-edge spar is spliced at the center line and has no cap. The 
ribs are segmented, being broken at the spars. The cover sheets are of 
uniform thickness and are continuous across the center line. The stringers 
are made of equal-legged angles, are continuous across the center line, 
and are riveted to the outside of the covers for convenience in construc- 
tion. The entire structure is symmetrical in the spanwise and depthwise 
directions and there are no cutouts . 
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CALCULATION OF GENERALIZED INFLUENCE COEFFICIENTS 


In the calculations, the approach outlined in reference 3 was 
followed in detail. The end product is a matrix which gives the deflec- 
tion (pQ, chordwise slope cp^, and chordwise half- cm* vat ure qp^ of 

neutral surface at a number of equally spaced stations along the trailing 
edge (y-axis) in terms of the total load Pq, moment p^, and second 

moment p^ at each station. This generalized influence -coefficient 

matrix [g] is obtained from the inverse of a matrix [a] which is 

formulated from the stiffness properties of the various parts of the 
structure. The calculations leading to the formulation of Qa] are 
illustrated herein by following the step-by-step procedure in the section 
entitled "Mechanics of Application" in reference 3* The numbers in paren- 
theses refer to the numbers of the steps in that paper, and the reader 
is advised to follow those steps in detail in conjunction with the present 
paper. All numerical values are in inch-pound units and are used with 
more significant figures than the actual measurements contain in order 
to avoid accumulated round-off errors. 

(l), (2) The chosen coordinate system and stations are indicated in 
figure 1. The y-axis lies along the center line of the web of the rear 
spar and the stations are spaced at seven l6-inch intervals. Thus the 
stations coincide very closely with the rivet lines of every second rib. 
Note that the positive direction of the x-axis is opposite to that of 
reference 3 * This change was made in order that points on the structure 
would have positive values of x. No change in notation is necessary 
except that c-^y) and C 2 (y) denote the x-position of the trailing 

edge and leading (swept) edge, respectively, instead of the reverse. Thus 

c x (y) = 0 

c 2 (y) =96 (0 < y < 16) 

c 2 (y) = 112 - y (16 < y < 112) 

(3) Because of symmetry, the neutral surface coincides with the 
middle surface. The stiffness properties of the spars and ribs are given 
in tables 2 and 3 of reference 1. The leading-edge spar (spar 5) must 
be treated as two separate spars 5(a) and 5(b)- The area of the stringers, 
the z-coordinate of the centroid of the stringers, the thickness of the 
cover sheets, and the z-coordinate of the midplane of the cover sheets 
are given in table 1 of reference 1. The x-coordinates of the stringers 
can be found from figure 2 of reference 1 if 5-25 inches is used as the 
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distance between the center line of the rear spar and the rivet line of 
the rearmost stringer. 

(4) The required values of the cover stiffness parameters n 

and a n are given in table I. They and all succeeding values of 

k,n+± 

2 

stiffness have been divided by E, which factor will ultimately reappear 
as a multiplier on the generalized influence-coefficient matrix. 


(5)* (6), (7) The values of effective stiffness P g n and location 

x s for the swept spar 5(h) are given in table II. The contributions 

of the unswept spars and the stringers are included in the combined 
fashion suggested in step (lO) of reference 5. Consequently, the summed 


quantities 


} s ,n* 


^s,n x s, n> 


, n x s. 


n 7 


I 


, n x s 


n 


and 


I 


3 


4 

s,n x s , n 


are listed in table II. 


(8) The values of the rib stiffness parameter 7 for each rib are 
given in table III; the ribs are numbered in accordance with their 
station locations, either on or halfway between stations. It should be 
noted that the rib stiffnesses have been reduced somewhat (about 5 P er ~ 
cent) below those in reference 1 in order to account for the ribs* being 
broken at the spars. 

From the foregoing quantities, the matrix [a] can be set up, 
inverted, and modified to obtain the [ g] matrix. Tables IV(a) and IV(b) 
show the results for the cases of symmetric deformation and antisymmetric 
deformation, respectively. In either case, 

I <p| = [s]|p| (!) 

For symmetric deformation | cp | is a 24-element column matrix: eight 

(Pq's, eight cp 1 's, and eight cp 2 's for stat i° ns 0 through 7. The column 

matrix |p| is similarly made up of eight p^'s, P-^' s > and- Pg ' s • F° r 

antisymmetric deformation | cp j and | p | are 21-element matrices, the 
quantities at the center line not appearing because they are always zero. 
Note that the boundary conditions specified in reference 3 require that 
cp and cp be zero in the symmetric case and that cp be zero 

0,0 1,0 °>i 

in the antisymmetric case. Space has been provided for these quantities 
in the matrix formulation in order to allow the introduction of rigid-body 
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notions necessary for the analysis of free-free vibrations. At this 
point, the rows of zeros in the [ g] matrices are sufficient to satisfy 
the boundary conditions on the supported structure. 


DETERMINATION OF GENERALIZED MASS MATRIX 


In this section, the calculation of the generalized mass matrix 
needed for the determination of vibration modes and frequencies is out- 
lined. Since the mechanics of the application of the Stein-Sanders 
approach to vibration problems has not heretofore been described, a 
detailed treatment is desirable here. 

As has been pointed out in reference 3, the generalized loads p^ n 
are given in terms of the potential energy of the applied loads II p by 




' k, n 




k,n 


( 2 ) 


where, for sinusoidal natural vibration. 



Note that w = \|/ Q + xty, + x 2 i|> 2 j where the i|r k ’s are the actual gener- 
alized deflections which include possible rigid-body motions. The use 

of t, in equation (2), rather than <p, as was employed in refer- 
k,n K,n 

ence 3, is valid here since the relation between the inertia loading 

and the vibration amplitude is independent of the particular boundary 

conditions on the structure. Carrying out the integration in the chord- 

wise direction yields 


■ - i I‘[ 


ra oV 


+ m. 


.( 2 Vh) + m 2 ( 2 ^0^2 + '•'l 2 ) + 2® 3 V2 + "Vs 


&y 


(M 


where 


nii. 


■l 


mx^dx 


( 5 ) 
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The functions m^ are not, in general, continuous. Such things 

as structural discontinuities and fuel cells will produce finite dis- 
continuities in m^; concentrated masses and ribs will produce impulse- 

type discontinuities. It is convenient to assume temporarily that there 
are no discontinuities, and then to correct for their effects subse- 
quently. Thus, for continuous m^, trapezoidal integration of equa- 
tion (4) gives 


n p - ' ^r(f “OjO^OjO + m O,l' l 'o,l + • • • + m o,N-l’*'o,N-l + 2 “o, 1^0,1 


2g 


m l, 0^0, 0^1,0 + 1*1,1 + • • • + ^H-iVn-I+I.N-I + + 


| m t,0' tl 2,0 + m 4,1^2,l + • * • + m 4,N-1^2,N-l + 2 


(6) 


where the equally spaced stations are numbered from 0 at the center 
line to N at the tip. 


Using equation (2) yields 
fr 


CD 

IT 


2 m 0,0 


ern 


0,1 


iu 


€ 

— m n ~ 

2 1,0 


em- 


1,1 


2 m 2,0 


em 2, 1 


em 0,N-l 


2 m 0,N 


r 0 


em 


1,N-1 


| m l,N 


em, 


2,N-1 


2 “2^K 


r\ 


(7) 
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and and p 2 can be expressed in a similar fashion. 

Thus, 



or 



where the are (N + l) X (N + l) matrices. The contribution of the 

continuous n% to these matrices is given by 



A finite discontinuity in m^ can be handled in several ways, but 
the simplest is merely to alter the value of at the station nearest 

the discontinuity. This can be done in precisely the same manner as is 
outlined for discontinuous a^ in step (k) of the step-by-step proce- 
dure in reference 3 ■ 

The energy due to a concentrated mass or a rib mass located at a 
distance ed to the right of station j is 


✓ 

n p = - - d) Vj + d Vj + i ] 2 + ^ip 1 - dH o,j + d % J+ i] p 1 - d ^i ;J + + 

- dH 0,j + d V J+l] [ (1 - d)+ 2, j + d +2,j + x] + ^P 1 - d) *l,J + d+ l, j+l] 2 + 

^P 1 - d) *l,j + d +l,J+l] P 1 ' d) *2,j + d +2,j + l] + %P X - d >*2,J + d ^2,j + l] 2 } 
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where the m^/s are integrals across the chord analogous to those 
defining Note that linear interpolation has been used to specify 

the values of between stations. Using equation (2) gives 




p 0,j = (1 ■ d) Vo,J + d(l - d ) m o' t 'o,j+i + (1 - d )* m l*l,j + d ( 


2- 


l,j+l 


(1 - d) 2 m^ 2>J + d(l - d)m 2 ^ 2) J+1 


y do) 


...2 f 


P 0,j+1 = T^ (l ■ d) ”oVj + d ”0+0, J+1 + d(l - d)m !+i,j + d "1+1, J+1 + 


d (l * d >”2+2,j + d2m 2^2,j 


j+1 


and so forth. Thus it can be seen that the_ effects of the concentration 
can be incorporated by including in each J^M^J matrix a block of terms 
like 

(l - d) 2 in k d(l - d)m k 

d(l - d)m k d% k 

The upper left-hand element is to be located in the ( j + l) row and 
column. All the concentrations are handled in this manner. 

For the delta-wing specimen under consideration the generalized 
mass matrix is given in table V. The computed numbers are based on the 
mass (weight) data given in tables 2 , 4 , and 5 and figure 4 of refer- 
ence 1. The generalized mass matrix also includes the masses of shaker 
armatures (2.0 pounds each) and pickups (0.7 pound each) located as 
shown in figure 10 of reference 1. 

In order to determine free-free modes and frequencies, the con- 
straining conditions used in obtaining the influence coefficients must 
be relaxed. The freeing procedure leading to the calculation of these 
vibration characteristics is discussed in the next section. 

FREE -FREE VIBRATION ANALYSIS 


Symmetric Modes 

For symmetric free-free vibration the actual generalized deflec- 
tion 4 n is related to the constrained generalized deflection cp n by 
the following matrix equation: 
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1*1 - 1*1 + ♦o,ol I ol + *l,o| I l| < U > 

where J Iq | and | I-jJ are (for the present structure) 24-element columns. 
The first eight elements of |Iq| are equal to 1 and all others are 0; 
the second eight elements of |lj| are 1 and the rest are 0. The term 
^0,0|1()| thus allows rigid-body translation; the term i|r-^ q|IjJ & BBows 
rigid-body pitching. 

Expressing | qp | in terms of |p | (eq. (l)) and then |p| in terms 
of |i|r| (eq. (9)) yields 

1*1 - ^[s]MI*l + * 0,0 N + *l,o| I l| (12) 


The unknowns i|)q q and q can be found from the conditions for self- 
equilibration of the inertia loading; that is, 

Nip'-l (i3) 

LhJM - 0 

v 

where the rows [^IqJ and j^I^J are the transposes of the corresponding 
columns. Multiplying equation (12) by [m] and then by [IqJ and 
gives , respectively, 

0 = M [g] DOUI + *06*0,0 + B oi+i,o 
0 = ^EiJ M fe] MEI + B io+o,o + B nh,o 

where the left-hand sides are 0 by virtue of equation (13) and 

Boo - EoJ M M 


(14) 




B 01 “ B 10 


- foj M l J il 

E n - [iij M |iil 


(15) 


Note that Bqo half the mass of the specimen, Bq-j_ is half the mass 
moment of the specimen about the trailing edge, and B-^ is half the 
mass moment of inertia of the specimen about the trailing edge. 
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Equations (l4) can be solved for \|> 0 Q and ^ Q . The result is 


v 0,0 


* 1,0 


-f — 211 — ^r( B 0i FaJ - B 11 Foj) M M [M] 1 1 
s( B 00 B n - B oi) V y 




s ( B 00 B 11 B 01 


( B oiFoJ - Wij)M[g]Mbl 


> ( 16 ) 


Substitution of 0 811(1 0 into e ^ ua ‘ bion (12) yields, 

finally, 

1*1 - ^([i] + ps] }[£!][«] 1+ 1 (17) 

where [ij is the identity matrix and 

[%] ( B Ol|Iol blj + B 01 |Il| foj - 

B 00 B 11 " B 01 ^ 

B OoPllfiJ - ^xFol N)M (18) 

Equation ( 17 ) is in a suitable form to be handled by standard 
iteration techniques in order to calculate the natural modes and 
frequencies . 


Antisymmetric Modes 

For antisymmetrical free-free vibration, the following relationship 
exists: 


Id - M + * 0 ,iM ( 19 ) 

where |r| is (for the present structure) a 21-element column with the 
first seven elements equal to the integers from 1 to 7 and “the other 
elements equal to 0, The term q|rj allows freedom of rigid -body 

rolling. 

Using equations (l) and ( 9 ) gives 

M = ^[g] [M] UI + V Q}1 \r\ 


( 20 ) 
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Note that in this case [m] is a matrix of twenty-first order obtained 
from the one in table V by deleting the first, ninth, and seventeenth 
rows and columns of that matrix. It should be recalled that this reduc- 
tion of order comes about because the generalized deflections at the 
center line (station zero) are necessarily zero. 

Self-equilibrium of the inertia rolling moments requires that 

L r J I P I = 0 (21) 


or 



0 = fid MM MM + “*0,1 

(22) 

where 

p - LrJ [M]|r| 

(23) 

is the mass moment 

of the half -span about the center line. 


Solving for 

\|tq 2 _ and substituting in equation ( 20 ) gives 



hi = J- [m + Pa]} [eJMM 

( 24 ) 

where 

rr1 MhJM 

L f aJ = " n 

(25) 


Equation (24) is in a form suitable for iteration. 


RESULTS AND DISCUSSION 


The first four modes and frequencies for free -free symmetric and 
antisymmetric vibration were calculated by matrix iteration from equa- 
tions ( 17 ) and (24). The results are given in tables VI and VII and in 
figure 2. In table VI are shown the frequencies obtained from the Stein - 
Sanders method, the frequencies measured experimentally (ref. l) , and 
the percentage difference. In table VII are shown the calculated mode 
shapes. These are given in terms of the values of at each station 

for each mode. The node-line patterns for the various modes are shown in 
figure 2. The node-line locations at each station were calculated from 
the equation 

+ xtyj + x^2 = 0 
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Also shown for comparison are the node- line patterns obtained experi- 
mentally (ref. l). 

As shown by figure 2, the node lines calculated by the Stein-Sanders 
method agree extremely well with those obtained experimentally. On the 
other hand, the errors in the calculated frequencies shown in table VI 
are large, especially for the higher modes. These errors are most prob- 
ably due to inaccuracies in the generalized influence coefficients; these 
inaccuracies, in turn, arise from basic shortcomings of the Stein-Sanders 
method when applied to this particular structure. The most serious of 
these are the lack of enough freedom in the chordwise shape of the deflec- 
tion and the neglect of the effects of transverse shear. The magnitude 
of these effects can be visualized from the comparison of theoretical and 
experimental static deflections under a uniform load shown in figure 3* 

Figure 3 ( a ) shows the deflections along the odd-numbered stations 
due to a uniform loading of 1 psi for the specimen mounted on the three- 
point support used in the experiment (ref. l) . Figure 3(b) shows the 
deflections for the specimen cantilevered at the root chord. The theo- 
retical deflections were calculated from generalized influence coeffi- 
cients modified for the different support conditions by the procedure 
outlined in appendix A of reference 3* The experimental deflections 
were calculated from the experimental influence coefficients tabulated 
in reference 1. 

Figure 3 shows large errors in the theoretical deflections, partic- 
ularly in the region close to the center line. Here the curvatures in 
the chordwise direction are large. Going from the three-point support 
to the cantilever support greatly reduces the magnitude of the curvature 
in the chordwise direction and improves the theoretical results, as can 
be seen from the solid curves in figure 4. Further improvement results 
from including approximately the effects of transverse shear, as illus- 
trated by the long -and -short -dashed curves. The transverse- shear cor- 
rections, which were obtained merely by adding the shear deflections of 
the structure treated as a beam to the solid-line deflections, bring 
the theoretical results into excellent agreement with the experimental 
results . 


CONCLUDING REMARKS 


The poor agreement between theoretical and experimental vibration 
frequencies indicates that the Stein-Sanders method is unsatisfactory 
for the analysis of the particular delta- wing specimen treated herein. 
To conclude that the method is unsatisfactory for more realistic struc- 
tures, however, would be erroneous since the present specimen had no 



extra stiffening in the chordwise direction such as would be afforded by 

!^ lag V n “ aCtual case * The static-deflection results imply that 
a significant improvement in the accuracy would result from such stif- 
ening. On the other hand, substantial errors due to tranverse-shear 

effects would still occur unless these effects were incorporated into 
the analysis. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., October 22, 1958. 
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TABLE IV.- GENERALIZED INFLUENCE- COEFFICIENT MATRIX 
(a) Symmetrical Deformation 


600 

«01 

602 




g 10 

g ll 

g 12 

’ g iJ = g ji' 



g 20 

g 21 

622 





0 


0 

0 

0 



0 

1.4605 

2.3371 

2.9633 

3.4715 

3.9320 

4.3812 


0 

2.3371 

6.2959 

9.0016 

11.1183 

12.9979 

14.8160 

g 00 = 10 2 

0 

2.9633 

9.0016 

16.5142 

22.0266 

26.7453 

31.2486 

0 

3.4715 

11.1183 

22.0266 

34.8672 

44.9403 

54.2410 


0 

3.9320 

12.9979 

26.7453 

44.9403 

66.4296 

84.4999 


0 

4.3812 

14 . 8160 

31.2486 

54.2410 

84.4999 

121.7644 


0 

4.8134 

16.5740 

35.6145 

63.2448 

101.7525 

154.7403 



0 

0 

0 

0 

0 

0 

0 


0 

-5.2356 

-7.3041 

-8.1018 

-8.3574 

-8.4477 

-8.6696 


0 

-7.3528 

-18.4467 

-22.7682 

-24.1506 

-24.5853 

-25.3639 

1n -4 

0 

-8.2481 

-23.0549 

-40.2180 

-45-6460 

-47.0877 

-48.7870 

g oi = 10 

0 

-8.6298 

-24.8620 

-46.4386 

-71.4180 

-77.3319 

-80.6397 


0 

-8.8514 

-25.7434 

-48.8833 

-79-1175 

-116.6625 

-125.35IO 


0 

-9.0593 

-26.4959 

-50.6333 

-82.85OO 

-126.8145 

-191.9687 


0 

-9.2374 

-27.1494 

-52.1841 

-86.0003 

-133.1657 

-205-4808 


0 

0 

0 

0 

0 

0 

0 


-4.5894 

-0.3323 

0.9345 

0.9774 

0.5543 

0.1588 

0.6541 


-10.4855 

-4.8743 

3.0949 

4.2470 

2.5854 

0.7063 

2.0006 

a - 1 O - 6 

-16.4006 

-10.5071 

-0.7644 

10.1746 

7.8673 

2.5061 

3.8709 

g 02 “ 10 

-21.7915 

-16.0319 

-6.5318 

5.2467 

19.6448 

9.0176 

6.3516 


-26.5575 

-21.0319 

-12.3566 

-2 .4834 

9.9552 

30.4915 

10 . 8050 


-30.9116 

-25.6066 

-17.7943 

-10.1823 

-3.3790 

5.2410 

48.4020 


-35.0030 

-29.9262 

-22.9798 

-17-5840 

-16.4639 

-23.7835 

-44.7327 


0 

0 

0 

0 

0 

0 

0 


0 

28.3524 

37.9775 

40.3652 

39.9856 

39-1644 

39-1042 


0 

37.9775 

95.1750 

111.8504 

112.2741 

109.1631 

108.4191 

„ in -6 

0 

40.3652 

111.8504 

204.3501 

218.2796 

210.5924 

206.4018 

fill = 10 

0 

39.9856 

112.2741 

218.2796 

372.4092 

372.8470 

356.8939 


0 

39.1644 

109.1631 

210.5924 

372 . 8470 

667.7383 

634.3781 


0 

39 • 1042 

108.4191 

206.4018 

356.8939 

634.3781 

1,327.3634 


0 

39.0262 

108.0850 

205.3515 

353.0658 

616.0558 

1,236.7512 



0 

0 

0 

0 

0 

0 

0 


25.8394 

-1.4518 

-9.8037 

-9.9815 

-6.7716 

-3.3508 

-2.7615 


53.9159 

18.0223 

-35.8236 

-43 .4052 

-30.6770 

-14.8455 

-9.1998 

812 - io ~ 8 

78.6843 

41.3218 

-24.2851 

-115.1064 

-98.2423 

-51.0132 

-23.9552 

97.5257 

61.3070 

-2.7311 

-102.7624 

-275.8926 

-179.4049 

-68.4267 


109-6814 

74.8382 

14.8797 

-74.3519 

-242.4892 

-666.7485 

-253.1988 


117.1688 

82.8651 

25.0295 

-55.9268 

-193.0163 

-525.2459 

-1,954.3729 


116.9185 

82.6621 

24.9975 

-55.3751 

-189.3031 

-498.6749 

-1,808.3827 


159.4650 

131.6396 

100.0772 

69.7104 

43.9338 

24.5883 

6.5715 


131.6396 

132.8083 

110.6374 

80.9280 

51-9840 

28.9754 

9.3431 


100.0772 

110.6374 

138.8273 

118.6105 

80.3315 

44.0332 

16.3945 

822 = 10“ 10 

69.7104 

80.9280 

118.6105 

210.1371 

171.7124 

97.7948 

37.5284 

43-9338 

51.9840 

80.3315 

171.7124 

459.4673 

324.9323 

122.8315 


24 . 5883 

28.9754 

44.0332 

97-7948 

324.9323 

1417.3669 

597.8610 


6.5715 

9.3431 

16.3945 

37-5284 

122.8315 

597.8610 

6,995.8497 


21.2506 

24.4307 

52.7155 

58.8470 

148.5812 

461.6689 

3,673.1493 


0 

4.8134 

16.5740 

33.6145 

63.2448 

101.7525 

154.7403 

230.3527 


0 

-8.6415 

-25.2526 

-48.4920 

-79.8611 

-122.8877 

-182.5040 

-286.4404 


0 

0.0363 

-0.4154 

-2.306l 

-9.0831 

-37.3008 

- 176.2185 

-733.9283 


0 

39.0262 

108.0850 

205.3515 

353.0658 

616.0558 

1,236.7512 

7,756.8221 


0 

-2.0778 

-6.4535 

-16.8594 

-40.1844 

-50.2044 

37.4970 

49,693.9410 


21.2506 

24.4307 

32.7155 

58.8470 

148.5812 

461.6689 

3,673.1493 

511,957.79 
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TABLE IV.- GENERALIZED INFLUENCE -COEFFICIENT MATRIX - Concluded 
(b) Antisymmetrical Deformation 


3 g 00 ®01 

[g] - y s 10 g n 

g 20 g 21 


0 0 0 

o 2.57125 4.36403 

0 4.56405 10.9165 

g = 10* 2 0 5.85200 15.9004 

0 7.20010 20,5374 

0 8.55431 24.6510 

0 9.91682 28.9694 


0 

0 

0 

1.79887 

-4.53310 

-6.08128 

4.61612 

-2.40976 

-14.3545 

8.03087 

2.74318 

- 11.0121 

11.7777 

9.04169 

- 3.65303 

15.6755 

15.6450 

4.56756 

19.6327 

22.3059 

12.8404 


g 02 

g 12 ; g ij = g Ji' 
g 22 


0 

0 

0 

0 

5.83200 

7.20010 

8.55431 

9.91682 

15.9004 

20.3374 

24.6510 

28. 9694 

28.5161 

38.7284 

48. 3060 

57.8162 

38.7284 

60.8982 

79.8337 

98.2626 

48 . 3060 

79.8337 

118.5971 

153.4663 

57.8162 

98.2626 

153.4663 

231.7661 

0 

0 

0 

0 


- 5.78413 - 5.15891 - 4.74374 - 4.72819 
-16.2288 -15.2881 - 14.4178 - 14.3140 
-30.4154 -32.4931 -31-6401 - 31.1675 
-26.2776 -58.5211 -61.5330 -59.4878 
- 17.9999 - 55.1910 - 113-1384 - 104.2097 
-9.42059 -48.4080 -112.0715 - 193.6732 


0 

0 

0 

0 

0 

0 

0 

0.0490400 

5.96449 

6.36279 

4.39499 

2.11223 

0.506115 

0.880214 

0.115012 

7.32222 

16.6890 

13.0266 

6.0558 

1.03498 

0.696020 

0.186129 

6.90149 

16.3851 

28 . 3729 

14.8358 

1.72992 

- 4.31509 

0.257123 

6.06107 

13.0697 

22.1589 

38.5729 

4.83003 

- 30.3993 

0.326777 

5.26017 

9.53703 

12.1339 

15.4922 

41.0693 

-167.0757 

0.395885 

4.53113 

6.22020 

2.38588 

-11.3255 

-53.2766 

-722.6711 


12.7550 
20.4739 
c 25.7624 

g ll = 10 29.5069 

32.0013 
33.6062 

33.4978 


20.4739 25.7624 
73.8791 91.9924 
91.9924 183.3723 
95.6147 198.1404 
95.2772 191.5197 
95.2198 186.0352 
94.9139 185.1218 


29.5069 32.0013 
95.6147 95.2772 
198.1404 191.5197 
350.7493 350.4543 
350.4543 643.2705 

330.9064 605.0374 
327.3572 587.0817 


33.6062 33.4978 

95.2198 9^.9139 

186.0352 185.1218 

330.9064 327.3572 

605.0374 587.0817 

1 , 291.820 1 , 201.730 

1,201.730 7,722.306 


0.357139 

0.416251 

0.372644 

g = 10'° 0.304273 

0.250619 

0.222539 

0.222330 


2.51667 4.23019 
- 35.7211 - 35.9270 
- 42.3980 -H9.9726 
- 37.4011 -H9.2887 
- 31.1834 - 98.9304 
- 27.4213 - 84.4719 
- 27.3427 - 83.9134 


4.88730 4.53501 

-22.0109 -7*57499 
- 96.6797 - 46.7637 
-281.2345 -I78.II68 
- 252.5471 -667.5713 
-204.4838 -526.4627 
- 200.8741 - 500.0124 


2.48129 6.66447 

-0.649586 3.97644 

-11.9076 -6.58982 

- 52.0014 -30.0105 

- 235.6573 - 41.7332 

- 1 , 935-246 43.6259 

-1,789.626 -49,687.76 


-0.576591 

- 0.437242 

-0.364781 

- 0.245316 

-O.136065 

-0.0606687 

- 0.103730 

-0.457212 

43.2112 

47.7417 

34.7671 

18. 8111 

6.86534 

11.6145 

- 0.364781 

47.7417 

157.5055 

138.0318 

79.2549 

28.5229 

45.0942 

- 0.245316 

34.7671 

138.0318 

438.0681 

313.1202 

113.7636 

137.7855 

-O.I36065 

18.8111 

79.2549 

313.1202 

1 , 411.054 

593.1698 

457.5185 

-O.O606687 

6 . 86534 

28 . 5229 

113.7636 

593.1698 

6,992.212 

3,676.822 

- 0.103730 

11.6145 

43.0942 

137.7855 

457.5185 

3,676.822 

511 , 958.0 



TABLE V.- GENERALIZED MASS MATRIX 


20 
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TABLE VI.- COMPARISON OF VIBRATION FREQUENCIES 


Mode 

Symmetric modes 


Antisymmetric modes 

Frequency, cps 

Error, 

Frequency, cps 

Error, 


Experimental 

Theoretical 

percent 

Experimental 

Theoretical 

percent 

First 

43.3 

46.4 

7.16 

52.2 

56.70 

8.43 

Second 

88.8 

105.3 

18.58 

91.7 

103.4 

12.77 

Third 

122.8 

150.0 

i 

22.15 

131.1 

166.6 

27.08 

Fourth 

164.2 

202.0 

23.02 

169.2 

216.5 

27.95 
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TABLE VII.- THE t COEFFICIENTS FOR SYMMETRIC AND ANTISYMMETRIC MODES 


Station, 

First symmetric mode 

First antisymmetric mode 

n 

CJ 

+l,n 

+2,1! 

+0,n 

+l,n 

+2,n 

0 

- 1.4862 

0.0214 57 

- 0.00005546 

0 

0 

0 

1 

-1.3789 

.023320 

-.00005404 

-.0312 

.005752 

.00000017 

2 

-1.0651 

.020134 

-.00005258 

-.1154 

.013131 

-.00001392 

3 

-.3229 

.014939 

- .00005700 

-.2687 

.021039 

- .00002562 

4 

.2582 

.007699 

-.00007280 

-.5137 

.028935 

-.00002058 

3 

1.2596 

-.001205 

-.00010434 

-.8625 

.036493 

.00001118 

6 

2 .4220 

-.011701 

- .00009519 

-1.2995 

.043408 

.00006029 

7 

3.6672 

-.014848 

-.00069212 

-1.7965 

.045575 

.00050698 


Station, 

Second symmetric mode 

Second antisymmetric mode 

n 

+0,n 

+l,n 

+2,n 

+0,n 

+l,n 

+2,n 

0 

0.2255 

- 0.009340 

0.00012448 

0 

0 

0 

1 

.2154 

-.010401 

.00013217 

.3181 

- .006477 

- .00000030 

2 

.1974 

-.013968 

.00015982 

.6018 

- .013888 

.00001913 

3 

.1804 

-.019120 

.00019049 

.7396 

- .019110 

.00004005 

4 

.1983 

-.025141 

. 00020134 

.6272 

- .019885 

.00005835 

5 

.2954 

-.031593 

.00016630 

.1908 

- .014828 

.00008543 

6 

.4930 

-.037493 

.00003626 

-.5527 

- .004515 

.00014013 

7 

.7923 

- .041011 

-.00052598 
i 

-1.5264 

.003366 

.00100000 


Station, 

n 

Third symmetric 

mode 

Third antisymmetric mode 

+0,n 

♦l,n 

+2,n 

+0,n 

+l,n 

♦a,n 

0 

-0.0682 

- 0.005360 

0.00009133 

0 

0 

0 

1 

.0204 

-.006366 

.00009484 

-.0643 

- .001502 

.00000042 

2 

.1029 

-.OO8858 

.00010193 

-.1578 

.002960 

-.00006281 

3 

.2486 

-.011115 

.OOOO9855 

-.2096 

.012590 

-.00018514 

4 

.3240 

-.011129 

.00007847 

- .2107 

.026305 

-.00033897 

5 

.2105 

-.006905 

.00005602 

-.2266 

.041626 

-.00043593 

6 

-.1625 

.001761 

.00009912 

-.3680 

.053467 

-.00021608 

7 

-.8011 

.010870 

.00093706 

-.7626 

.062547 

.00100000 


Station, 

n 

Fourth symmetric mode 

Fourth antisymmetric mode 

+ 0 ,n 

+l,n 

+ 2 ,n 

+ 0 ,n 

♦l,n 

+ 2 ,n 

0 

0.4506 

- 0.025574 

0.00025835 

0 

0 

0 

1 

.3874 

- .023308 

. 00023509 

-.0964 

.002070 

-.00000012 

2 

.1998 

-.014407 

.00015347 

- .1062 

.000798 

.00002517 

3 

-.0619 

.000218 

.00000601 

- .0097 

-.003736 

.00006835 

4 

-.2976 

.018320 

- .00018550 

.1522 

-.009523 

.00010543 

5 

-.3649 

.033068 

-.00033811 

.2373 

-.011913 

.00010629 

6 

-.1756 

.035897 

-.00028145 

.0574 

- .005948 

.00009922 

7 

.2303 

.029159 

- .00035149 

-.5128 

.007359 

.00096160 












3rd Mode 


4th Mode 


(b) Antisymmetric modes. 


Figure 2.- Node-line patterns for free- free natural modes. 
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